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ABSTRACT 

Context. NGC7538 is a complex massive star-forming region. The region is composed of several radio continuum sources, one of 
which is IRS 1, a high-mass protostar, from which a 0.3 pc molecular bipolar outflow was detected. Several maser species have been 
detected around IRS 1. The CH3OH masers have been suggested to trace a Keplerian-disk, while the H2O masers are almost aligned 
to the outflow. More recent results suggested that the region hosts a torus and potentially a disk, but with a different inclination than 
the Keplerian-disk that is supposed to be traced by the CH3OH masers. 

Aims. Tracing the magnetic field close to protostars is fundamental for determining the orientation of the disk/torus. Recent studies 
showed that during the protostellar phase of high-mass star formation the magnetic field is oriented along the outflows and around or 
on the surfaces of the disk/torus. The observations of polarized maser emissions at milliarcsecond resolution can make a crucial con- 
tribution to understanding the orientation of the magnetic field and, consequently, the orientation of the disk/torus in NGC7538-IRS 1. 
Methods. The NRAO Very Long Baseline Array was used to measure the linear polarization and the Zeeman-splitting of the 22 GHz 
H2O masers toward NGC7538-IRS 1. The European VLBI Network and the MERLIN telescopes were used to measure the linear 
polarization and the Zeeman-splitting of the 6.7 GHz CH3OH masers toward the same region. 

Results. We detected 17 H2O masers and 49 CH3OH masers at high angular resolution. We detected linear polarization emission to- 
ward two H2O masers and toward twenty CH3OH masers. The CH3OH masers, most of which only show a core structure, seem to trace 
rotating and potentially infalling gas in the inner part of a torus. Significant Zeeman-splitting was measured in three CH3OH masers. 
No significant (3(t) magnetic field strength was measured using the H2O masers. We also propose a new description of the structure 
of the NGC7538-IRS 1 maser region. 

Key words. Stars: formation - masers: water, methanol - polarization - magnetic fields - ISM: individual: NGC7538 



J^ [ 1 . Introduction pling of the gas with a large-scale magnetic field around IRS 1. 

■ . Because Sandell et al. (2009) found that the collimated free-free 

I^ NGC7538 is a complex massive star-forming region located j^j (opening angle < 30°) is approximately aligned with the out- 

O in the Perseus arm of our Galaxy at a distance of 2.65 kpc fl^^ ^nd that there is a strong accretion flow toward IRS 1 (ac- 

^H (Moscadelh et al. 2009 ). The region is composed of seve ral clus- ^retion rate ~ 2 x 10^ Moyr'), IRS 1 must be surrounded by 

^ , ters of infrared sources (Wynn- William s et al. 11974)1 and ra- ^n accretion disk. The morphology of the free-free emission, 

•^ ■ dio continuum sources (CampbeU[198l. The brightest source ^^ich is opticall y thick up to 100 GHz (Franco-Hernandez & 

. ^ IS NGC7538-IRS 1, whose central star has been suggested to Rodriguez |200l, suggests that the disk should be almost edge- 

^ be an 06 star of about 30 Mq with systemic local stand ard of on and oriented eas t-west (ScoviUe et al 2QQ8; Kameya et al 

H rest velocity Visr = - 58 kms' (Campbell & Thompson [19841 [19891 SandeU et al. [20091 . A possible detection of this edge- 

. .^ , Sandell et al. '2009', Puga et al. ,2010,). Several high-velocity on disk was made by Minier et al. ( 1998"), who observed a linear 

molecular bipolar outflows were detected in NGC7538 and one distribution of the brightest CH3OH masers with an inclination 

of these IS elongated 0.3 pc from IRS 1 (position angle PA=140°) ^ngle of about 1 12°. Pestalozzi et al. (l2004l) estimated that this 

with a veloc ity of 250 km s' aii d a mass of 8 2.8 M p (Kameya jisk is a Keplerian disk with an outer radius of -750 AU and an 

et al [19891 Gaume et al. [1995] Davis et al. [1998[ Qm e t al. ^^^^j. radius of -290 AU by modeling the CH3OH maser emis- 

1201 11 1. VLA continuum observations by Campbell ( |1984| l in- sjons at 6 7 and 12 2-GHz 

dicate that the PA of the outflow decreases away from IRS 1; Recent results disagree with the edge-on disk ti-aced by the 
i.e., 1 80° at ^.'3 (0.004 pc) and 165° at 2" (0.03 pc). Kameya CH3OH masers. The observations of the mid-infrai-ed emis- 
et al. ([198911 gave three possible interpretations of this rotation: ^jon suggest that the radio continuum emission does not come 
disk precession, interaction of flie flow with dense gas, and cou- f^om a free-free jet but that it traces the ionized gas wind from 
the disk surface that in the new scenario would be perpendic- 



* Member of the International Max Planck Research School (IMPRS) ular to the CO-bipolar outflow with a disk inclination angle 
for Astronomy and Astrophysics at the Universities of Bonn and / = 32° (De Buizer & Minier [2005l l. Klaassen et al. (120091) ob- 
Cologne. 
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Fig. 1. Positions of water, methanol, and hydroxyl maser features superimposed on the 2 cm continuum contour map of NGC7538- 
IRS 1 observed with the VLA in 1986 (Franco-Hernandez & Rodriguez I2004I I . Contours are 1, 2, 4, 8, 16 xl.3 mJybeam '. The 
(0,0) position is aaooo = 23'^13'"45'*.382 and ^2000 - +61°28'10'.'441. Blue circles and red triangles indicate the positions of the 
H2O and CH3OH maser features, respectively, with their linear polarization vectors (60 mas correspond to a linear polarization 
fraction of 1%). Green boxes indicate the positions of the OH masers detected by Hutawarakorn & Cohen (2003) with their linear 
polarization vectors (thick lines, 6 mas correspond to 1%). The two arrows indicate the direction of the bipolar outflow (~140°). 



served two warm gas tracers (SO2 and OCS) toward NGC7538- 
IRS 1 with the Submillimeter AiTay (SMA). Although the re- 
gion was unresolved, they found a velocity gradient consistent 
with the CH3OH maser velocities and perpendicular to the large- 
scale molecular bipolar outflow. This rotating gas might indicate 
that there is a torus (with an angular size of about 2 arcsec; i.e., 
~5300 AU at 2.65 kpc) surrounding the smaller accretion disk 
proposed by De Buizer & Minier (2005]l. 

In addition, Krauss et al. (I2006I I proposed two possible sce- 
narios for NGC7538-IRS 1 by considering both the linear dis- 
tribution of the 6. 7 -GHz CH3OH masers detected previously 
(e.g., Pestalozzi et al. 2004 De Buizer & Minier 12005b and 
the asymmetry in the near infrared (NIR) images observed 
by them. In their most likely scenario, called "Scenario B", 



the CH3OH masers trace the edge-on disk as suggested by 
Pestalozzi et al. (I2004I I and the asymmetry might be caused by 
the precession of tiie jet. In the "Scenario A" the CH3OH masers 
do not trace a disk but an outflow cavity as proposed by De 
Buizer & Minier ( 20051 1. where the detected asymmetry might 
simply reflect the innermost walls of this cavity. 

Besides the CH3OH masers, other maser species were de- 
tected around NGC7538-IR S 1: O H, H2O, NH3, and H2CO 
(e.g., Hutawa rakorn & Cohen 12003] Galva n-Madrid et al. l20T0l 
Gaume et al. 1991', Hoffman et al. I2003I I. The OH masers are 
located southward and show no obvious disk structure or re- 
lation to the outflow direction (Hutawarakorn & Cohen 2003JI. 
H2CO and H2O masers are located near the center of the contin- 
uum emission, and the H2O masers are also almost aligned with 
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Table 1. All 22 GHz H2O maser features detected in NGC7538-IRS 1. 
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Notes. ^"^ The best-fitting results obtained by using a model based on the radiative transfer theory of HtO masers for F + F^ - Is"' (Surcis et al. 12011 1 . The errors were determined by 
analyzing the full probability distribution function. "^' The percentage of circular polarization is given by Py = (Vmax ~ Vniin)/^iax- ''^'The angle between the magnetic field and the maser 
propagation direction is determined by using the observed P\ and the fitted emerging brightness temperature. The errors were determined by analyzing the full probability distribution 
function. 



the outflow (Galvan-Madrid et al. l20T0l l. The 6.7 and 12.2-GHz 
CH3OH masers show a cone shape that opens to the north-west 
(Minier et al. 2000). 

So far, the magnetic field structure in NGC7538-IRS 1 has 
been studied using submillimeter imaging polarimetry (Momose 
et al. I2OOII 1 and OH maser emission (Hutawarakorn & Cohen 



12003 ). Although the polarization vectors are locally disturbed, 
at an angular resolution of 14 arcsec the magnetic field direc- 
tions agree with the direction of the outflow (Momose et al. 
[2001), while at milliarcsecond (mas) resolution the OH maser 
observations indicate a magnetic field oriented orthogonal to 
the outflow (Hutawarakorn & Cohen 120031 ) . Because H2O and 
CH3OH masers were detected close to the center of the contin- 
uum emission, i.e. to the protostar, it is worthwhile investigat- 
ing their linear and circular polarization emissions. As shown by 
Surcis et al. (2009 2011) for the massive star-forming region 
W75N, the polarization observations of the two maser species 
can depict a reasonable scenario for the magnetic field close to 
massive protostars. Moreover, the direction and strength of the 
magnetic field might help to decide the debate about the orienta- 
tion of the disk. 

Here we present Very Long Baseline Array (VLBA) 
observations of H2O masers, Multi-Element Radio Linked 
Interferometer network (MERLIN) and European VLBI 
Network (EVN) observations of CH3OH masers in full po- 
larization toward NGC7538-IRS 1. In Sect. |3] we show our 
results obtained by studying the linear and circular polariza- 
tion of H2O and CH3OH maser emissions in a similar way as 
Vlemmings et al. (2006, 2010) and Surcis et al. (2009, 201 1 ) for 
Cepheus A and W75N, respectively. In Sect. |4] we discuss our 
results and attempt to disentangle the complex morphology in 
NGC7538-IRS1. 



2. Observations and analysis 

2.1. 22 GHz VLBA data 

We observed the star-forming region NGC7538-IRS1 in the 616- 
523 transition of H2O (rest frequency: 22.23508 GHz) with the 



NRACQ VLBA on November 21, 2005. The observations were 
made in full polarization spectral mode using four overlapped 
baseband filters of 1 MHz to cover a total velocity range of « 50 
kms '. Two correlations were performed. One with 128 chan- 
nels to generate all four polarization combinations (RR, LL, RL, 
LR) with a channel width of 7.8 kHz (0.1 kms '). The other 
one with high spectral resolution (512 channels; 1 .96 kHz=0.027 
kms '), which only contains the circular polarization combi- 
nations (LL, RR), to be able to detect Zeeman-splitting of the 
H2O maser across the entire velocity range. Including the over- 
heads, the total observation time was 8 h. 

The data were edited and calibrated using the Astronomical 
Image Processing System (AIPS) following the method of 
Kemball et al. (I1995I ). The bandpass, the delay, the phase, and 
the polarization calibration were performed on the calibrator 
J0359-(-5057. The fringe-fitting and the self-calibration were per- 
formed on the brightest maser feature (W06 in Table [1]). All 
calibration steps were initially performed on the dataset with 
modest spectral resolution after which the solutions, with the 
exception of the bandpass solutions that were obtained sepa- 
rately, were copied and applied to the high spectral resolution 
dataset. Stokes /, Q, and U data cubes (4 arcsec x 4 arcsec, 
rms ^ 7 mJy beam"') were created using the AIPS task IMAGR 
(beam-size 1 .0 mas x 0.4 mas) from the modest spectral reso- 
lution dataset, while the / and V cubes (rms ^ 8 mJybeam"') 
where imaged from the high spectral resolution dataset and for 
the same fields. The Q and U cubes were combined to produce 
cubes of polarized intensity and polarization angle. Because 
these observations were obtained between two VLA polarization 
calibration observation^ made by the NRAO in 2009, during 
which the linear polarization angle of J0359+5057 was constant 
at -86?7, we were able to estimate the polarization angles with a 
systemic error of no more than ~3°. 

We identified the H2O maser features using the process de- 
scribed in Surcis et al. (i201II ). Here the program called "maser 
finder" is used. This program is able to search maser spots, ve- 



' The National Radio Astronomy Observatory (NRAO) is a facility of 
the National Science Foundation operated under cooperative agreement 
by Associated Universities, Inc. 

" http://www.aoc.nrao.edu/ smyers/calibration/ 
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Fig. 2. Left panel: a close-up view of the H2O maser features around NGC7538-IRS1. Right panel: a zoom-in view of the boxed 
region of the left panel. The octagonal symbols are the identified maser features in the present work scaled logarithmically according 
to their peak flux density. The maser LSR radial velocity is indicated by color. A 10 Jybeam ' symbol is plotted for illustration 
in both panels. The reference position is aaooo - 23*'13™45^.350 and ^2000 - +61°28'10'.'428, which was estimated as described in 
Sect. [3] The linear polarization vectors, scaled logarithmically according to polarization fraction Pi (in Table [T]|, are overplotted. 
The dashed line (PA=-52°) is the result of the best linear fit of the H2O masers of group N (features W16 and W17 were not 
included). The synthesized beam is 1.0 mas x 0.4 mas. 



locity channel by velocity channel. We identified a maser fea- 
ture when three or more maser spots coincide spatially (within 
a box 2 by 2 pixels) and each of them appeared in consecutive 
velocity channels. In Table [T] only the brightest spot of each se- 
ries of maser spots that fulfill the criteria described above are 
reported. The two maser features that show hnear polarization 
emission were fitted using a full radiative transfer method code 
based on the models for H2O masers of Nedoluha & Watson 
(II992I 1. They solved the transfer equations for the polarized ra- 
diation of 22 GHz H2O masers in the presence of a magnetic 
field which causes a Zeeman-splitting that is much smaller than 
spectral line breadth. The fit provides the emerging brightness 
temperature (T\,AQ.) and the intrinsic thermal linewidth (AVi). 
See Vlemmings et al. C006b and Surcis et al. (1201 11 1 for more 
details. We modeled the observed linear polarized and total in- 
tensity maser spectra by gridding AVi from 0.5 to 3.5 kms"', 
in steps of 0.025 kms"', by using a least-squares fitting routine. 
From the fit results we were able to determine the best values of 
the angle between the maser propagation direction and the mag- 
netic field {0). By considering the values of TbAQ, AVi, and we 
were also able to estimate the saturation state of the H2O maser 
features. 



2.2. 6. 7 GHz MERLIN data 

To detect the polarization of the CH3OH maser emission at 
6668.518MHz (5i - 60 A+) we observed NGC7538-IRS1 with 
six of the MERLIl^ telescopes in full polarization spectral mode 
on December 28, 2005. The observation time was ~ 2 h, includ- 
ing overheads on the calibrators 2300-1-638, 3C84 and 3C286. 
We used a 250kHz bandwidth (~ llkms ') with 256 channels 
(velocity resolution ~ 0.04 kms"') centered on the source ve- 



locity Visr = -56.1 kms '. For calibration purposes, the con- 
tinuum calibrators were observed with the 16 MHz wide-band 
mode. Both 3C84 and 3C286 were also observed in the narrow- 
band spectral line configuration and were used to determine the 
flux and bandpass calibration solutions. The data were edited 
and calibrated using AIPS. The calibrator 3C84 was used to de- 
termine the phase off'set between the wide- and narrow -band set- 
up. Instrumental feed polarization was determined using the un- 
polarized calibrator 3C84, and the polarization angle was cal- 
ibrated using 3C286. Using one of the strongest isolated maser 
features, we were able to self-calibrate the data in right- and left- 
circular polarization separately. After calibration, the antenna 
contributions were reweighed according to their sensitivity at 
5 GHz and their individual efficiency. Stokes /, Q, U data cubes 
(5.12 arcsec x 5.12 arcsec, rms s; 10 mJybeam"') were created 
(beam-size 47 mas x 34 mas). 

The CH3OH maser features were also identified through the 
identification process mentioned above (Surcis et al. 1201 11 1. If 
the program finds a group of maser spots that are not exactly 
spatially coincident but show a continuum linear distribution 
with a clear velocity gradient, we report in the corresponding 
table only the brightest maser spot of the group (e.g. maser fea- 
ture M06). We were unable to identify weak CH3OH maser fea- 
tures (< 1 Jybeam"') close to the brightest ones because of 
the dynamic range limits. To determine TbAQ, AVi, and 6, we 
adapted the code used for 22 GHz H2O masers to model the 
6.7 GHz CH3OH masers, which has successfully been used by 
Vlemmings et al. (2010) for the CH3OH masers in Cepheus A. 
We modeled the observed linear polarized and total intensity 
CH3OH maser feature spectra by gridding AVi from 0.5 to 
2.6 kms"', in steps of 0.05 kms"', by using a least-squares fit- 
ting routine. 



' MERLIN is operated by the University of Manchester as a National 
Facility of the Science and Technology Facilities Council 
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Table 2. All 6.7 GHz CH^OH maser features detected in NGC7538-IRS 1 with MERLIN. 
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0.86 
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+39 ±8 


< 0.55" 


> 12" 


- 


M02 


A 


-115 


47 


1.70 ±0.02 


-55.94 


0.45 


- 


- 


- 


- 


- 


M03 
M04 


A 
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0.63 
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61.83 ±0.01 
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0.4 ±0.1 
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+ 18±13 
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Notes. '"^^ The absolute reference position is 0-2000 - 23'M 3"" 45^.362 and ^2000 - +61°28'10'!506. *'^* The best-fitting results obtained by using a model based on the radiative transfer theory 
of CH3OH masers for F + F,, = ls~' (Vlemmings et al. fTOiO i The errors were determined by analyzing the full probability distribution function. '''The angle between the magnetic field 
and the maser propagation direction is determined by using the observed P\ and the fitted emerging brightness temperature. The eiTors were determined by analyzing the full probability 
distribution function. '^'Because of the low angular resolution the code is able to give an upper limit for AVj and a lower limit for TbAQ. 



2.3. 6.7 GHz EVN data 

NGC7538-IRS1 was also observed at 6.7 GHz (CH3OH) in 
full polarization spectral mode with nine of the EV>0 anten- 
nas (Jodrell2, Cambridge, Efelsberg, Onsala, Medicina, Torun, 
Noto, Westerbork, and the new joint antenna Yebes-40 m), for 
a total observation time of 5 h on November 3, 2009 (program 
code ES063B). The bandwidth was 2 MHz, providing a veloc- 
ity range of ~ 100 kms"'. The data were correlated using 1024 
channels to generate all four polarization combinations (RR, LL, 
RL, LR) with a spectral resolution of 1.9 kHz (~0.1 kms '). 

The data were edited and calibrated using AIPS. The band- 
pass, the delay, the phase, and the polarization calibration were 
performed on the calibrator 3C286. Fringe-fitting and self- 
calibration were performed on the brightest maser feature (E26 
in Table O. Then we imaged the /, Q, U, RR, and LL cubes 
(2 arcsec X 2 arcsec, rms a; 8 mJybeam"') using the AIPS task 
IMAGR (beam-size 6.3 mas x 4.9 mas). The Q and U cubes 
were combined to produce cubes of polarized intensity and po- 
larization angle. 

In this case, because the dynamic range of our EVN obser- 
vations was better than that of MERLIN observations, we were 
able to identify CH3OH maser features with a peak flux den- 
sity of less than 1 Jybeam"' with the "maser finder" program. 
The maser emission was fitted by using the adapted code for 
CH3OH masers, but with a grid of AVi ranging from 0.5 to 
1.95 kms '. We were able to determine the Zeeman-splitting 
from the cross-correlation between the RR and LL spectra, 
which was successfully used by Surcis et al. (2009) for the po- 
larized CH3OH maser emission detected in W75N. The dynamic 
range of the RR and LL cubes decreases close to the strongest 
maser emission of each group because of the residual calibration 
errors. As a result, we were not able to determine the Zeeman- 
splitting (AVz) for the features with a peak flux density of less 
than 1.8 Jy beam '. 



3. Results 

In Fig. [1] we show the H2O (blue circles) and CH3OH maser 
features (red triangles) detected by us and superimposed on 



'' The European VLBI Network is a joint facility of European, 
Chinese, South African and other radio astronomy institutes funded by 
their national research councils. 



the 2 cm continuum contour map of NGC7538-IRS1 observed 
with the VLA in 1986 (Franco-Hernandez & Rodriguez 2004). 
Because the continuum observations were made 20.3 yr be- 
fore our observations, we have shifted the continuum map by 
-50 mas in both directions according to the proper motion, 
fia - -2.45 mas yr"' and //^ = -2.45 mas yr"', measured by 
Moscadelli et al. (120091 ). The direction of the large-scale molec- 
ular bipolar outflow is also shown. Because we did not have ab- 
solute positions of the H2O maser features, we estimated the 
offset of one common H2O maser feature detected both in the 
VLA observations of Galvan-Madrid et al. (2010) and by us (fea- 
tures M3 and WOl, respectively). All positions of the H2O maser 
features were shifted according to this offset. We did indeed 
have absolute positions of the CH3OH maser features at low- 
angular resolution, but not at high-angular resolution. Hence, 
to overlay these maser features to the continuum, we matched 
the brightest CH3OH maser feature detected with the EVN and 
the brightest one detected with MERLIN (M06, absolute posi- 
tion 0-2000 = 23''13'M5^.362 and 52000 = +6r28'10':506). AU 
CH3OH maser features were shifted accordingly. Thus, based 
on the resolution of the VLA H2O and MERLIN CH3OH maser 
observations the uncertainties of H2O and CH3OH maser fea- 
tures absolute positions are 50 mas (Galvan-Madrid et al. 1201 Oi l 
and 10 mas, respectively. 



3.1. H2O masers 

We detected 17 22-GHz H2O maser features with the VLBA as- 
sociated to NGC7538-IRS1 (named W01-W17 in Table [B. No 
H2O maser emission with a peak flux density (/) less than 0.13 
Jybeam"' is detected even though our channel rms is signifi- 
cantly less. 

The H2O maser features can be divided into two groups, N 
and S, that are composed of 1 1 and 6 H2O maser features, re- 
spectively. Group N is located at the center of the continuum 
emission and it shows a linear distribution with a position angle 
of -52° (see Fig.|2]i. The local standard-of -rest velocities (Vki) of 
group N are between -60.2 km s"' and -43.3 km s"' . Excluding 
W16 and W17, the range is -60.2 kms-i< V!'* < -58.3 kms"'. 

' '=' isr 

Group S is located about 1" southward from group N and its 
velocity range is -73.3 kms"'< V^^^ < -57.3 kms"', which is 
more blue-shifted. 

Linear polarization is detected in 2 H2O maser features that 
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Fig. 4. Left panel: CH3OH maser structures detected with MERLIN around NGC7538-IRSL Contours are 16, 32, 64, 128, 256, 
512, 1024 X 0.15 Jybeam"'. The synthesized beam is 47 mas x 34 mas. Right panel: a close-up view of the CH3OH maser 
features detected with the EVN around NGC7538-IRS1. The synthesized beam is 6.3 mas x 4.9 mas. The triangle symbols are 
the identified maser features in the present work scaled logarithmically according to their peak flux density. The maser LSR radial 
velocity is indicated by color. A 5 Jy beam"' symbol is plotted for illustration in the right panels. The absolute reference position is 
aim.) - 23*^1 3'"451362 and ^aooo = +61°28'10'.'506 as determined from the MERLIN observations. The linear polarization vectors, 
scaled logarithmically according to the polarization fraction Pi (Tables |2] and |3]l, are overplotted. 




AVJkms'] 

Fig. 3. Results of the full radiative transfer ;^^^-model fits for 
the H2O maser feature W07. The fit yields the emerging maser 
brightness temperature TbAQ and the intrinsic maser ther- 
mal linewidth AVi. Contours indicate the significance intervals 
h.x^-0-25, 0.5, 1, 2, 3, 7, with the thick solid contours indicating 
Icr and 3cr areas (see VIemmings et al. 120061 and Surcis et al. 
I2011l for more details). 



belong to group S (W06 and W07). Their linear polarization 
fraction (column 7 of Table [TJ is Pi » 1 % and the weighted 
linear polarization angles is (x^ q) ~ -11°. The full radiative 
transfer method code for H2O masers was able to fit the feature 
W07 and the results are given in column 10 and 11 of Table [T] 
and in Fig. [3] The emerging brightness temperature and the in- 
trinsic thermal linewidth are lO'' K sr and 3.4 kms ', respec- 
tively. By considering T\,AQ. and the observed Pi we determined 
1°, indicating that the maser is operating in a regime 



gB 

'H2O 



where the magnetic field is close to perpendicular to the propa- 



gation of the maser radiation. No significant circular polarization 
emission is detected. 



3.2. CH3OH masers 

We detected 13 6.7-GHz CH3OH maser features with the 
MERLIN telescope (named MO 1 -Ml 3 in Table |2]l that appear 
to be composed of 49 features when observed with the EVN 
resolution (named E01-E49 in Table [3]). Including the more sen- 
sitive EVN observations, we detect no CH3OH maser emission 
below 0.1 Jy beam"'. Note that as defined in Sect. [2] we indicate 
with the term maser feature the brightest maser spot among a 
series of maser spots that either show a spatial coincidence and 
consecutive velocities or show a clear velocity gradient along a 
continuum linear structure. 

In the left panel of Fig. |4] we show the contours of the 
CH3OH maser structures detected with MERLIN and in the right 
panel the distribution of the CH3OH maser features detected 
with the EVN. The CH3OH maser features distribution at high 
angular resolution match the CH3OH maser emission detected 
with MERLIN four years before perfectly. Following the nam- 
ing convention adopted by Minier et al. (2000), they can be di- 
vided into five groups (from A to E). Note that each group is 
composed of several maser features, each of which indicates a 
series of maser spots. Group A, which is composed of five maser 
features (3 at MERLIN resolution) and hosts the brightest maser 
feature of the region (M06 and E26, respectively), shows a linear 
distribution. If we consider all maser spots of the maser feature 
M06 (top-left panel of Fig. |5]l, we are able to measure a velocity 
gradient of about 0.02 kms"' AU"' from northwest to southeast, 
which is confirmed by considering the matching maser spots of 
E26 (EVN, top-right panel of Fig.©. 

Group B is resolved into 12 CH3OH maser features (four 
at MERLIN resolution). The spots of the features E33 and E40 
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Fig. 5. Zoom-in view of the CH3OH maser spots of the features M06, which was detected with MERLIN, E26, E33, and E40, which 
were detected with the EVN. Features M06 and E26 are the same CH3OH maser detected with the two instruments. The LSR radial 
velocity of the spots is indicated by color Ten, 4, and 2 Jy beam ' symbols are plotted for illustration in the panels. 



of this group show a velocity gradient similar to that of E26 
but from southwest to northeast (bottom panels of Fig. |5]l. No 
CH3OH maser features of group B detected with MERLIN show 
a velocity gradient. The CH3OH maser groups C, D, and E are 
composed of 18, 3, and 11 maser features, respectively. While 
groups C and D are located close to the central peak of the 
continuum emission, group E is located about 300 mas south- 
ward. The velocities of groups C and D are more blue-shifted 
(-61.5 kms"' < Visi < -60.5 kms"') than those of groups A, B, 
andE(-59kms-i < Vi^ < -56 kms"'). 

We detected linear polarization in 10 and 20 CH3OH maser 
features with MERLIN and EVN, respectively. The features of 
group D show the highest linear polarization fraction of the re- 
gion (column 8 of Tables |2] and [3), though the observations at 
higher angular resolution revealed E43 to be the feature with 
the highest linear polarization fraction (6.2%). Because the EVN 
provides an angular resolution eight times better than MERLIN 
and because we also have a higher dynamic range, we will, for 
the interpretation of the magnetic field, only use the linear po- 
larization vectors of the maser features detected with the EVN. 
The groups A, C, and E have weighted linear polarization vec- 
tors almost oriented east-west, with angles -86°+31°, -88°±6°, 
and -i-89°±44°, respectively. The weighted linear polarization 
angles for the other two groups are (Xqh ob) ~ -60°±36° and 

(^?H30h) - +48°±32°. 

The full radiative transfer method code for CH3OH masers 
was able to fit 7 (MERLIN) and 20 (EVN) CH3OH maser fea- 




Fig. 6. Result of the full radiative transfer ;if^-model fits for the 
CH3OH maser E33 detected with the EVN. The fit yields the 
emerging maser brightness temperature T\,AQ. and the intrinsic 
maser thermal linewidth AVi. Contours indicate the significance 
intervals A;t'^=0.25, 0.5, 1, 2, 3, 7, with the thick solid contours 
indicating Icr and 3cr areas. 



tures, the results are given in columns 10 and 11 of Tables |2] 
and |3] Considering the high angular resolution observations, 
the weighted intrinsic maser linewidth and the weighted emerg- 
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Table 3. All 6.7 GHz CH^OH maser features detected in NGC7538-IRS 1 with the EVN. 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(H) 


(12) (13) 


Ma,ser 


group 


RA 

offset 


Dec 

offset 


Peak flux 
Density(I) 


Vkr 


AVL 


Pt 


X 


Ay° 


7-1, An" 


AV| 9'' 






(mas) 


(mas) 


(Jybeam"') 


(kms"') 


(kms"') 


(%) 


(°) 


(kms"') 


(log K sr) 


(ms-') (■■) 


EOl 


D 


-286.910 


-282.502 


0.55 ±0.02 


-60.40 


0.12 


- 


- 


- 


- 


- 


E02 

E03'' 

E04 


D 
D 
C 


-272.725 
-255.344 
-93.660 


-279.037 
-273.992 
-222.357 


16.82 ±0.03 
4.72 ± 0.03 
10.39 ± 0.02 


-60.49 
-60.58 
-61.28 


0.21 
0.24 
0.24 


4.5 ±0.8 

5.6 ±0.9 
3.6 ±0.2 


+ 13 ±16 

+58 ±5 
+81 ±22 


1.0+°-' 
0,9+0:5 

LO+I 


94J+0.94 

'■-"'-0.47 

9 7fi+°-5* 
-0.14 


+2.7 + 0.3 77+] J 
80+?„ 
83+?3 


EOS 


C 


-86.171 


-200.790 


0.73 ± 0.02 


-61.28 


0.20 


- 


- 






- 


E06 
E07 
E08 
E09 
ElO 


A 
C 
A 
C 
C 


-80.377 
-74.392 
-61.547 
-57.174 
-52.938 


22.786 
-213.730 

18.163 
-207.658 
-248.180 


15.21 ±0.03 
16.41 ±0.03 
4.08 ± 0.02 
9.48 ± 0.02 
0.73 ± 0.02 


-57.24 
-61.28 
-57.15 
-61.46 
-60.84 


0.20 
0.29 
0.23 
0.27 
0.31 


3.6 ±0.2 
4.1 ±0.5 
4.3 ±1.4 
2.5 ±0.7 


-66 ± 10 
-83 ± 20 
-48 ± 40 
-90 ± 4 


1 0+»-2 
0.8!| 


Q -74+0.69 

939?8:il 
9 3Y+o:6i 
9.10+*f 


77i 
79+^1 


Ell 


c 


-49.085 


-205.700 


1.12 ±0.03 


-60.75 


1.08 


- 


- 


- 


- 


- 


E12 


c 


-48.483 


-245.093 


1.21 ±0.02 


-60.67 


0.28 


- 


- 


- 


- 


- 


E13 
E14 


c 
c 


-44.930 
-41.350 


-222.357 
-198.801 


4.65 ± 0.02 
1.11 ±0.01 


-60.84 
-61.02 


0.25 
0.49 


1.4 ±0.4 


-86 ± 83 


1 2+°-' 

' -^-0.3 


8.8 r°S 


'°-38 


EI5 


c 


-40.421 


-255.891 


0.59 ± 0.02 


-60.84 


0.33 


- 


- 


- 


- 


- 


E16 


c 


-39.082 


-181.780 


0.57 ± 0.03 


-60.75 


0.18 


- 


- 


- 


- 


- 


EI7 
EI8 


A 

c 


-37.661 
-27.986 


12.650 
-199.362 


14.31 ±0.02 
0.80 ±0.03 


-56.80 
-60.75 


0.33 
0.19 


1.4 ±0.2 


+58 ± 12 


1 6+°-' 

1 ■°-0.4 


8.77^3 


8o:iJ! 


E19 


c 


-25.963 


-230.701 


1.29 ±0.02 


-60.84 


0.26 


- 


- 


- 


- 


- 


E20 
E21 


c 
c 


-22.083 
-18.147 


-211.354 
-226.401 


3. 19 ±0.03 
2.53 ±0.03 


-60.58 
-60.75 


0.27 
0.22 


2.4 ±1.0 
2.41 ±0.04 


-87 ± 23 
-83 ±12 


1 0+°' 


9.05+°^^ 
9 ns+o:56 

-0.06 


79+j; 

85!^; 


E22 


c 


-15.769 


-242.950 


0.51 ±0.02 


-60.84 


0.25 


- 


- 








E23 


A 


-13.091 


-56.523 


1.43 ±0.05 


-55.75 


0.26 


- 


- 


- 


- 


- 


E24 


c 


-8.172 


-220.826 


1.83 ±0.02 


-60.67 


0.29 


2.5 ±0.3 


-84 ±14 


l-3!Si 


9m*"„f, 


82-39 


E25 


c 


-7.762 


-238.071 


0.42 ± 0.02 


-60.84 


0.32 


- 


- 




- 




E26'' 


A 








95. 15 ±0.08 


-55.92 


0.33 


5.83 ±0.03 


-87 ±17 


0.5+°;* 


11 39+°" 

-0.03 


+ 1.6 + 0.3 86+5, 


E27 


E 


71.113 


-503.039 


0.19 ±0.01 


-58.56 


0.12 


- 


- 






- 


E28 


E 


81.252 


-497.234 


0.37 ±0.01 


-58.56 


0.27 


- 


- 


- 


- 


- 


E29 


E 


85.297 


-557.914 


1.06 ±0.01 


-58.91 


0.27 


- 


- 


- 


- 


- 


E30 


E 


97.923 


-513.329 


0.30 ±0.01 


-58.38 


0.18 


- 


- 


- 


- 


- 


E31 


E 


99.399 


-489.943 


2.26 ±0.01 


-58.38 


0.21 


- 


- 


- 


- 


- 


E32 
E33 


B 
B 


110.386 
116.016 


-235.856 
-175.760 


4.37 ± 0.02 
10.11 ±0.02 


-57.77 
-57.86 


0.31 
0.30 


2.4 ±0.6 

2.5 ±0.7 


+74 ± 44 
-81 ±83 


1 4+°-' 

-0.4 


9.06+?*; 
9.05+°« 


89+!° 
76!« 


E34 


B 


121.782 


-203.613 


0.54 ± 0.02 


-57.42 


0.55 


- 


- 




- 




E35 


E 


122.083 


-511.884 


0.80 ± 0.02 


-58.21 


0.22 


- 


- 


- 


- 


- 


E36 


B 


123.859 


-173.165 


6.46 ± 0.02 


-57.59 


0.32 


2.3 ±0.6 


-28 ± 42 


i-5!S:J 


9-04!|;™ 


89!]^ 


E37 


B 


125.062 


-233.976 


0.68 ± 0.02 


-57.07 


0.28 


- 


- 








E38 


E 


126.893 


-540.110 


0.93 ± 0.03 


-58.12 


0.21 


- 


- 


- 


- 


- 


E39 


E 


127.303 


-546.631 


0.21 ±0.06 


-58.73 


0.26 


- 


- 


_ 


- 


- 


E40 


B 


129.927 


-218.961 


8.98 ± 0.02 


-57.51 


0.37 


2.1 ±0.7 


-69 ± 10 


l-8-ol 


8.99+]°° 


^i-l 


E41'' 
E42 


B 

E 


132.550 
134.190 


-163.392 
-506.261 


1.81 ±0.02 
1.10 ±0.03 


-57.24 
-58.12 


0.34 
0.20 


4.3 ±1.1 


+76 ± 42 


o.6!s;= 


9.44+°™ 


85!5, 


E43'' 


B 


137.196 


-209.274 


1.69 ±0.03 


-57.24 


0.37 


6.2 ±2.0 


-18±17 


0.5+°;^ 


9.83+°^2 


90+™ 


E44 


B 


137.552 


-181.682 


1.05 ±0.02 


-57.33 


0.24 


- 


- 




- 


- 


E45 


B 


143.509 


-164.852 


0.98 ± 0.02 


-56.72 


0.32 


- 


- 


- 


- 


- 


E46 


E 


144.083 


-537.047 


0.12 ±0.01 


-58.73 


0.28 


- 


- 


- 


- 


- 


E47 


B 


149.112 


-203.584 


0.71 ±0.02 


-57.42 


0.27 


- 


- 


- 


- 


- 


E48 


E 


156.163 


-503.223 


23.68 ± 0.03 


-58.03 


0.24 


3.5 ±0.6 


+89 ± 44 


1.2!»;5 


9.26!i;»? 


-2.7 + 0.3 75+;^ 


E49 


B 


162.558 


-154.095 


1.50 ±0.05 


-56.19 


0.16 


- 


- 






- 



Notes. ''^' The best-fitting results obtained by using a model based on the radiative transfer theory of CH3OH masers for F + Fi, = Is"^ (Vlemmings et al. 2010 1. The errors were determined 
by analyzing the full probability distribution function. For F + F^ - 0.6s"^ (Minier et al. i2002t T\,AQ. has to be adjusted by adding -0.22. '''' The Zeeman-splittings are determined from 
the cross-conelation between the RR and LL spectra. *'^'The angle between the magnetic field and the maser propagation direction is determined by using the observed P[ and the fitted 
emerging brightness temperature. The enors were determined by analyzing the full probability distribution function. ''''Because of the degree of the saturation of these H2O masers, T^AQ. 
is underestimated, AV; and 9 are overestimated. 

Table 4. Weighted values of the linear polarization angles, the intrinsic thermal linewidths, the emerging brightness temperatures, 
and the angles between the line of sight and the magnetic field of the CH3OH masers for feature in common. 



Group 
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EVN 


(T,A(l) 


(0) 




MERLIN 






(AVi) 


(x) 


(AVi) 


(Tt,ACl) 


<e> 







(kms-') 


(Ksr) 








(kms-') 


(K sr) 


(°) 


A 


-86 + 31 


0.9!" 
0.5! 
1.0! 
0.6! 


4 
i 


5 X 10'-' 


86!J, 


-51+62 


<1.2° 


> 5 x lO"" 


90-53 


B 


-60 ± 36 


1x10' 


87+* 
85!^ 
79!f, 


-63 ± 76 


<0.9'' 


> 1 X lO'^* 




C 

D 


-88 ±6 

+48 ± 32 


i 


2x10' 
2x10'° 


-74 ± 69 
+39 ±8 


2 2+°' 
<0.6° 


3x 10" 
> 1 X 10'^* 


^61ll 


E 


+89 + 44 


1 2+»2 

1-^-0 4 


2x10' 


75!- 


+21 ± 54 


1 5+°> 


6x 10* 


m% 



Notes. *"* Only the highest value has been taken into account. ''^' Only the lowest value has been taken into account. 



ing brightness temperature are {Ayi)cH30H = 0.9 kms'' and 
(rbAQ)cH30H ~ 10^ K sr. As an example, the ;^f^-contours for 
feature E33 are reported In Fig. |6l 



The fit for the CH3OH maser features detected with 
MERLIN gives values both for AVi and for T\,AQ. consistent or 
higher than those detected with the EVN (see Table|4|i. However, 
owing to the lower angular resolution of MERLIN and the occur- 
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rence of strong velocity gradients, the fits are strongly affected 
by line blending. As for the H2O maser features, we are able 
to determine the 6 values (column 12 and 13 of Tables |2] and 
[3l respectively) from T\,ls.Q. and Pi. The weighted value for the 
whole region is (0ch,oh)evn - 86°;^43., which is almost constant 
in all groups. Only for group E we determine a lower value, i.e., 
<^ch,oh)evn =75°+^^„. 

The Zeeman-splitting (AVz) in ms ' determined from 
the cross-correlation between the RR and LL spectra of the 
CH3OH maser features detected with the EVN is reported in 
column 12 of Table [3] Note that the cross-correlation method 
is dynamic -range sensitive, therefore we are able to measure 
AVz only for those features that show a very high dynamic 
range. It has been impossible to obtain AVz measurements from 
MERLIN data because of the insufficient dynamic range. 

4. Discussion 

4.1. Comparing CH^OH spectra at different resolution 



Vlemmings ( I2008I I observed the CH3OH masers with the 100-m 
Effelsberg telescope in November 2007, i.e. two years later than 
our MERLIN observations and two years before our EVN ob- 
servations. A comparison between the single-dish and MERLIN 
fluxes reveals only a difference of < 10% across the whole spec- 
tra. This slight difference can be explained by the different res- 
olution of the two instruments and flux calibration uncertain- 
ties in the Effelsberg observations (~ 10%, Vlemmings 120081 ). 
This means that no significant changes occurred between the 
two epochs and there is no indication of a core/halo structure 
of the maser features (hereafter simply called masers) at scales 
of ~100 AU such as described for a sample of CH3OH maser 
sources by Pandian et al. ( 120^1 Ij i. 

In Fig. |7] we show the total flux spectra of the 




V|„ [km/sj 

Fig. 7. Profile of the total flux spectra of the CH3OH masers in 
NGC7538 detected with MERLIN (solid line, epoch 2005) and 
the EVN (dashed line, epoch 2009). The groups are also reported 
according to the velocity of Table [3] 

CH3OH masers detected with MERLIN and the EVN in 2005 
and 2009, respectively. Generally, the EVN resolved out be- 
tween 60% and 90% of the flux when 6.7-GHz CH3OH masers 



were observed with both these instruments (e.g., Minier et al. 
I2OO2 I Bartkiewicz et al.'2009'; Torstensson et aL'SoTl). This in- 
dicates the presence of emission structures on scales that are re- 
solved by the EVN, i.e. a halo structure that surrounds the core 
of the masers. Moreover, if the maser spots do not have a com- 
pact core, they would not be detected in the EVN observations 
(Pandian et al. 201 1 ). Although in our observations the EVN re- 
solves out about 50% of the flux of group C (Visr > -61kms"'), 
the majority of the flux appears to be recovered across the entire 
CH3OH maser spectrum. This suggests that the halo structure 
is absent, except in some masers of group C and likely in E26 
of group A, and that the masers only contain a compact core. 
Moreover, the non-detection of the maser M07 with the EVN 
might indicate the absence of the compact core in this maser. 

Minier et al. ('2002') observed the 6.7-GHz CH3OH masers 
with five antennas of the EVN. The authors determined the size 
of the halo (di) and of the core (d2) of the brightest masers 
of groups A and C by Gaussian fitting of the visibility ampli- 
tudes. For the brightest maser of group A the authors determined 
di - 5 mas (~13 AU at 2.65 kpc) and d2 - 3 mas (~8 AU), and 
di = 17 mas (~45 AU), and d2 = 5 mas (~13 AU) were mea- 
sured for the brightest maser of group C. The Gaussian fit of 
E26 and E07 gives a size of 5 mas and 6 mas respectively, indi- 
cating that the EVN resolves out the halo structure observed by 
Minier et al. (I2002I I and that we observe only the core of these 
two masers. 

We suggest that the absence of a halo structure in most of 
the CH3OH masers could be caused by the amplification of 
the strong continuum emission shown in Fig. [1] This would ar- 
gue that almost all CH3OH masers are on the front-side of this 
source. Similar arguments were used in the interpretation of the 
structure of Cepheus A (Torstensson et aI. I20III ). 

4.2. H2O and CH3OH maser properties 

Before discussing the polarization analysis of the masers in 
NGC7538-IRSI, we need to consider the degree of the maser 
saturation. As explained in detail by Vlemmings et al. (120 1 01) 
and Surcis et al. (201 1 ), when the full radiative transfer method 
code is applied to a saturated maser, the code gives a lower limit 
for TbAO and an upper limit for AVi. The masers are unsaturated 
whenR/(r + Ty) < 1 and fully saturated when /^/(r-Hrv) ~ 100, 
where R is the stimulated emission rate given by 



R 



AkBT^AO. 

Anhv 



(1) 



Here A is the Einstein coefficient for the maser transition, which 
for H2O masers is taken to be equal to 2 x 10^ s ' (Goldreich & 
Keelev [19721 and f or the CH3OH masers is 0.1532 x 10** s' 
(Vlemmings et al. 120101 ). ^b and h are the Boltzmann and 
Planck constants, respectively, and v the maser frequency. From 
RI(X + Ty) < 1 we can estimate an upper limit for the emerging 
brightness temperature below which the masers can be consid- 
ered unsaturated. The limits are (rbAQ)H,o < 6.7 x 10^ K sr and 
(rbAfi)cH30H < 2.6 X lO"* K sr for H26 and CH3OH masers, 
respectively. Consequently, only four CH3OH masers are partly 
saturated; i.e., E03, E26, E41 and E43 (Table [2). Because the 
model emerging brightness temperature scales linearly with F + 
Ty, the ratio RI(T + Ty) is independent of the value of F -H Fy. The 
Avl of the saturated CH3OH masers are close to AVi implying 
that the maser lines are rebroadened as expected when the maser 
becomes saturated. Furthermore, these four CH3OH masers do 
also show a high linear polarization fraction, which again con- 
firms their saturated state (Goldreich et aI. II973l) . 
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The intrinsic linewidth of the H2O maser W07 is higher 
than the typical values measured in previous works, i.e. AVi = 
2.5 km s"' (Surcis et al. 2011). This difference might indicate the 
presence of a strong turbulent gas with a turbulence velocity of 
AVtmb - 2.3 kms ' or multiple components overlapping. 

By comparing the brightness temperature T\, with TbAQ ob- 
tained from the model, we can estimate the maser beaming an- 
gle (AQ) for both maser species. We can estimate the brightness 
temperature (T\,) by considering the equation 



[K] 



S(y) 



[Jy] \ [mas^] 



■^, 



(2) 



where S{v) is the flux density, 2 the maser angular size and ^ 
is a constant factor that includes all constant values, such as the 
Boltzmann constant, the wavelength, and the proportionality fac- 
tor obtained for a Gaussian shape by Burns et al. (il979i ). The 
values of ^ for H2O and CH3OH masers are 



^22 GHz 



1.24- lO'^mas^Jy"' K, 



^6.7 GHz = 13.63-10^ mas^ Jy"' K, 



(3) 
(4) 



as reported by Surcis et al. (1201 11 1 and Surcis et al. ( I2009I I, re- 
spectively. The H2O maser W07 is unresolved and AQhjO - 
10 ~. The Gaussian fit of CH3OH masers detected with the 
EVN gives a size between 5 and 7 mas, which indicates that 
all CH3OH masers are marginally resolved. 

In a tubular geometry AD. » (d/l)^, where d and / are the 
transverse size and length of the tube, respectively. By assuming 
d approximately the size of the masers, the maser lengths are in 
the range lO'-' cm < 1^20 < 10'^ cm and 10'"^ cm < IchjOu < 
10'^ cm for W07 and the CH3OH masers, respectively. 

4.3. Magnetic field in NGC7538-IRS 1 

4.3.1. CH3OH Zeeman-splitting and magnetic field 

We measured 6.7-GHz CH3OH maser Zeeman-splitting ranging 
from -2.7 ms"' to -1-2.7 ms~^ This is 2-3 times larger than was 
measured with the Effelsberg telescope (Vlemmings I2008l l as 
expected when resolving individual masers. From the Zeeman 
theory we know that the Zeeman-splitting is related to the mag- 
netic field strength along the line of sight (By) by the equation 



AVz = ffz ■ Bn, 



(5) 



where AVz is the Zeeman-splitting in kms ' and az is the 
Zeeman-splitting coefficient in kms ' G~', which strongly de- 
pends on the Lande ^-factor of the molecular transition. 
Recently, Vlemmings et al. (1201 11 1 found an unfortunate cal- 
culation error in the az used so far (az = 0.049 kms ' G '; 
e.g., Surcis et al. '2009' Vlemmings et al. '2010"), the current 
best value is found to be one order of magnitude lower (az - 
0.005 kms"' G"'). However, Vlemmings et al. (2011) show that 
the splitting between the RCP- and LCP-signal still likely orig- 
inates from the Zeeman theory and not from other non-Zeeman 
effects. Until careful laboratory measurements of the ^-factor ap- 
propriate for the 6. 7 -GHz CH3OH masers are made, we cannot 
give any exact value for the magnetic field strength. Because 
there is a linear proportionality between AVz and B|(, we can 
say that the negative value of the Zeeman-splitting indicates a 
magnetic field oriented toward the observer and positive away 
from the observer. However, we can speculatively give a possi- 
ble range of values for |B||| by considering 0.005 kms"' G"'< 
az < 0.049 kms-' G"', thus 50 mG<|B||| < 500 mG. 



4.3.2. Magnetic field orientation 

Before discussing the orientation of the magnetic field, we have 
to evaluate the foreground Faraday rotation (i.e., the rotation 
caused by the medium between the source and the observer), 
which is given by 



Off] =4.22x10*' 



[kpc]j i[i 



Bn 



cm-3]/ \[mG]/ \[GHz] 



-2 



,(6) 



where D is the length of the path over which the Faraday rotation 
occurs, rtf. and By are respectively the average electron density 
and the magnetic field along this path and v is the frequency. By 
assuming the interstellar electron density, the magnetic field, and 
the distance are n^ * 0.012 cm-^ By » 2yuG (Sun et al. 2008l l. 
and D = 2.65 kpc, respectively. Of is estimated to be 0?5 at 22- 
GHz and 6?0 at 6.7-GHz. At both frequencies the foreground 
Faraday rotation is within the errors reported in Tables [1] and |3] 
thus it should not affect our conclusions. 

Because all 6 values are higher than the Van Vleck angle 
6'crit ~ 55° (Goldreich et al. I1973I I. the magnetic fields are in- 
ferred to be perpendicular to the linear polarization vectors. 
The magnetic field orientation derived from the H2O masers is 
^H,o = -1-79°, while the five groups of CH3OH masers show ori- 
entation angles <f^^^Q^ = +4°, ip^^^^^^ = +30°, <fi^^^^^^ = +2°, 

V'cH OH ~ -42°, and v:'ch,oh - ~^°- However, considering the 
relatively large errors in column 1 3 of Table [3] we cannot rule 
out that the actual 6 values are below 55°, in which case the 
magnetic field is parallel to the linear polarization vectors. 

Because Hutawarakorn & Cohen (120031 1 measured the lin- 
ear polarized emission of the 1.6 and 1.7-GHz OH masers, we 
can now compare the orientation of the magnetic field derived 
from the two maser species (Fig. [TJ. We have to consider that 
the Faraday rotation at the OH maser frequencies is, at the same 
conditions, larger than at 6.7-GHz. The linear polarization an- 
gles of group D show a 90° difference w.rt those of the nearby 
1.6-GHz OH masers. This might be caused either by a 90°-flip 
phenomenon or by the Faraday rotation if the masers are deeply 
located in a strong ionized gas. Most of the OH masers are 
likely to be Zeeman cr-components, therefore the magnetic field 
is perpendicular to the linear polarization vectors. Considering 
the large errors of 0e48> the orientation of the magnetic fields of 
group B and E are consistent with those derived from the nearby 
1.6-GHz OH masers (Hutawarakorn & Cohen 2003). In particu- 
lar, the linear polarization vector of E43, which shows the high- 
est linear polarization fraction, is perfectly aligned with feature 
2 of the 1665 MHz OH maser (Hutawarakorn & Cohen l2003l l. 
Unless the Faraday rotation is 180°, the consistence of the mag- 
netic field orientation indicates that the Faraday rotation at OH 
maser frequency is low. 



4.3.3. The role of the magnetic field 

The importance of the magnetic field in the region can be es- 
timated by evaluating the ratio between thermal and magnetic 
energies (J3). If p < 1, the magnetic field dominates the ener- 
gies in the high-density protostellar environment. The yS factor is 
given by 



/? = 2 



nif, 



(7) 
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where OTj is the Alfvenic Mach number and Ws is the sonic Mach 
number, which in formula are 



V3 



Va 



and 



m. 



V3 



Cs 



(8) 



Here Va is the Alfven velocity and Cs is the sound velocity, 
and o", the turbulence velocity, can be estimated using cr - 

AVtuib/ V8 ln2. By considering a mass /vh, = 3.8 x 10"^^ g and 
the relation between the velocity and the temperature of a gas, 
we can write Va and c^ in terms of |B|, tin,, and the kinetic tem- 
perature of the gas T^ 



Va 



[kms-i] 



[kms-i] 
So we get 



1542 



"H, 



\B\ 
[mG]/ \[cm-3] 



0.0603 



[K] 



/? = 3.058 -10 



-9 



\B\ Y^ I «H, 



[mG]/ \[cm-3]/ \[K] 



(9) 



(10) 



(11) 



which, considering Eq. |5] and \B\ = \B\\\/cos {0}, can thus be 
written as 



/3 ^ 3.058 ■ IQ-'^ al cos (0) 



"H; 



[ms-i]/ \[cm-3]/ \[K] 



3,,,^l-(12) 



By assuming mht = 10 cm , T ~ 200 K, which are the typical 
values in the CH3OH masing region, the weighted values for 
the unsaturated masers of the 6 angle {0} =78°, and Zeeman- 
splittings lAVzl = 2.7 ms"\ we obtain 



/?=17.4-4. 



(13) 



Although the value of az is still uncertain, we can expect 
that this must be neither higher than the old value (az - 
0.049 kms 'G ', Vlemmings 2008) nor lower than the new 
value (az = 0.005 kms-'G"\ Vlemmings et al. ISOTTT) . 
Therefore, it is reasonable that j6 is between 4 ■ 10"** and 4 • 10"^. 
Consequently the magnetic field is dynamically important in this 
massive star-forming region. 

4.4. Structure of NGC7538-IRS 1 

Recently, large-scale elliptical configurations have been de- 
tected in significant CH3OH masers around high-mass proto- 
stars. These masers are thought to trace molecular rings (e.g., 
Bartkiewicz et al. I2009I I. However, most of these rings do not 
show signs of rotation, but instead the radial motions dominate, 
indicating that the masers are instead tracing infalling gas in the 
interface between the torus and the flow, e.g. Cepheus A HW2 
(Torstensson et al. 2011; Vlemmings et al. 120101) . In this light 
it is legitimate to reconsider whether the masers in NGC7538- 
IRS 1 could trace similar gas. 

In this work we have observed velocity gradients in three 
CH3OH masers, i.e., E26, £33, and E40. Whereas the veloc- 
ity gradient of E26 has already been observed by other authors 
(e.g., Minier et al. 119981 Pestalozzi et al. 2004), those of masers 
E33 and E40 have not been reported so far (Fig.|5]). All velocity 
gradients of E26, E33, and E40 are equal to 0.02 kms"' AU"'. 
Considering the orientation of the linear distribution of the maser 
spots of E26, E33, and E40, we find that all of them seem to point 



toward a common center. This suggests that E26, E33, and E40 
are likely tracing a gas with a radial motion, probably consistent 
with the interface of the torus and infalling gas. Consequently 
the maser spots of E26 are likely not tracing a Keplerian disk 
as proposed by Pestalozzi et al. (.200 4). The CH3OH masers of 
groups A, B, C, and D show a cone-shape distribution that opens 
to the north-west. This seems to indicate that we are looking at a 
regular structure, for instance the inner part of a torus. Klaassen 
et al. ( 120091 ) detected an almost face-on 5300 AU-torus perpen- 
dicular to the outflow with a clockwise rotation motion. If this 
torus is on the same plane of the disk supposed by De Buizer 
& Minier (2005), then its inclination angle must be i=32°. By 
keeping in mind these new results, we can suppose that the 
CH3OH masers of groups A, B, and C are tracing the surface 
of the inner torus, possibly where the infall reaches the rotating 
structure. The masers of groups D and E trace the gas slightly 
farther out on the torus. In particular, group A is located to the 
opposite side of the torus w.r.t. the other three maser groups. 
In Fig. [8] a three-dimensional sketch of the region as described 
above is shown. The outflow sketched here is a simplified rep- 
resentation of the chaotic and complex multi-outflows structure 
observed toward NGC7538 and in particular around IRS-1 (e.g., 
Qiu et al. 201 1 , Klaassen et al. l20TTl ). 

In this scenario the H2O masers of group N are located par- 
allel to the edges of the outflow, which show a blue-shifted part 
northwestern and a red-shifted part southeastern (e.g., Kameya 
et al. 119891 Hutawarakorn & Cohen 120031 ). In this case the 
H2O masers can be pumped by a shock caused by the interac- 
tion of the outflow with the infalling gas. The most southern 
H2O masers, W06 and W07, might be associated either with the 
red-shifted part of the outflow or most likely with another source, 
as also suggested by the diff'erent orientation of the magnetic 
field in that region. Moreover, the OH masers located southward 
(Hutawarakorn & Cohen 120031 ) might be pumped by the red- 
shifted part of the high- velocity outflow and those located west- 
ward by the blue-shifted part. 

This scenario is further supported by the similarity between 
the torus and outflow velocities and that of the three maser 
species. In Fig.|9]we plot the CH3OH maser velocities as a func- 
tion of the position angle (PA) along the torus. To determine 
the PA of each individual CH3OH maser, we consider that ev- 
ery maser lies on an own circle that has the same orientation and 
inclination angle of the torus with the center located at the origin 
of the sketched outflow in Fig [8] Moreover, we plot in Fig. |9]an 
empirical function (solid line) based on the orientation and on 
the velocity field of the torus observed by Klaassen et al. ( 120091 ). 
i.e. the PA and velocity of all the points of a circle that has the 
same size, orientation, and inclination of the torus. The maser 
velocities perfectly match the non-Keplerian velocity profile of 
the torus that is rotating clockwise (Klaassen et al. 2009). This 
indicates that the CH3OH masers are likely related to the torus 
structure. Groups A and E appear to be systematically blue- and 
red-shifted w.r.t the torus velocities respectively, supporting the 
suggestion of infall motion (on the order of 1 kms"') along the 
line of sight. 

The velocities of the blue- and red-shifted part of the outflow 

VJ^d^ = -37 kms"' (Kameya et al. 

1989). The mean velocities of the H2O masers of group N and S 

59.4 kms"', excluding the masers W16 and W17, 

' . While y^ Q again suggests that these 
masers might be associated with another source, the difference 
between V^ q and V^^^^ indicates that the H2O masers are trac- 
ing the gas surrounding and being entrained by the outflow. The 



areV^i;; = -76 kms 

are V^ ^ 

and V^ o ~ ~68.7 kms 
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Fig. 8. The three-dimensional sketch of the massive star-forming region NGC7538-IRS 1 as described in Sect. 14.41 The two cones 
are the red- and blue-shifted part of the large-scale molecular bipolar outflow (Kameya et al. 1989 ), the donut is the torus suggested 
by Klaassen et al. ( I2009I I perpendicular to the outflow and with i=32° (De Buizer & Minier 2005 ), the triangles, the circles, and 
the boxes are the methanol, H2O masers (present paper), and the OH masers (Hutawarakorn and Cohen I2003I I. respectively. The 
outflow and the torus are transparent to allow the reader to see behind them. 



1.6 and 1.7-GHz OH masers show velocity ranges -61 kms"'< 
VoH < -51 kms"' that agree well with Vtoms and VchjOh, which 
might indicate an association with the torus rather than with 
the CO-outflow. The presence of OH and CH3OH masers in 
the same environments, also suggested by the linear polarization 
vectors of the two maser species (see Sect. 14.3.2b . is in contrast 
with the excitation model of Cragg et al. ( 2002) , which predicts 
the inhibition of the 6.7-GHz CH3OH masers when 1 .6-GHz OH 
masers arise. 

The orientation of the magnetic field derived by the linear 
polarized emission of the CH3OH masers at high angular res- 
olution matches the scenario we presented above. Most likely 
the magnetic field is perpendicular to the linear polarization vec- 
tors for most of the CH3OH maser features, but for group E the 
magnetic field is instead parallel to the linear polarization vec- 
tors (see Sect. 14.3.11 ). This is also suggested by considering the 
linear polarization vectors of the OH masers (see Sect. 14.3.21 ). 
In this case the magnetic field is situated on the two surfaces of 
the torus, with a counterclockwise direction on the top surface. 
This interpretation is consistent with the signs of the magnetic 
field strengths determined from the Zeeman-splitting measure- 



ments of the CH3OH masers and OH masers, which show neg- 
ative strength (B = -2.0 mG, Hutawarakorn & Cohen 120031) 
toward the region were groups B and E arises and positive to- 
ward groups A and D. 

Even if the CH3OH masers are not tracing the Keplerian-disk 
reported by Pestalozzi et al. (120041) . the presence of a smaller 
disk is not ruled out. In the scenario proposed here, it would be 
located within the torus. Accordingly, we argue that a possible 
answer to the question about the morphology of NGC7538-1RS 1 
is that the CH3OH masers are tracing a gas close to the torus that 
is falling toward the connection region between the torus and the 
disk as was observed in other protostars (e.g., Vlemmings et al. 
I2OTOI) . 

The new scenario described in this paper is determined by 
considering the entire CH3OH maser region and not just the 
masers of group A. Consequently, this scenario is different 
both from "Scenario A" and from "Scenario B" suggested by 
Kraus et al. 



(|2006)l, who considered the results obtained by De 
Buizer & Minier (2005) and by Pestalozzi et al. (2004) for the 
CH3OH masers of group A, respectively. In "Scenario A" the 
CH3OH masers of group A might resemble either clumps in the 
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Fig. 9. Measured velocities of the CH3OH masers as a function 
of the position angle (PA) along the torus. The solid line shows 
the velocities of the torus reported by Klaassen et al. (2009) as 
function of their PA along the torus determined as described in 
Sect. 14:41 



cavity of the outflow or recent ejecta from the outflow, while in 
"Scenario B" group A trace a disk and the authors proposed a 
jet precession model to explain the asymmetry of the NIR emis- 
sion. However, in our scenario the asymmetry of the NIR emis- 
sion might reflect the innermost walls of an outflow cavity as 
explained in "Scenario A", but with the masers that are tracing 
an infalling gas rather than the outflow cavity or a recent ejecta 
from the outflow. 

Two questions remain: what is the radio continuum emis- 
sion shown in Fig.[T]? How can this scenario be consistent with 
the changing of the position angle of the outflow reported in 
Campbell ( 1984)? De Buizer & Minier (I2005I I argue that the ra- 
dio continuum emission is consistent with radio emission aris- 
ing from a photoevaporated disk wind and Kameya et al. d 19891 1 
suggested that the rotation of the outflow stems from a disk pre- 
cession. 



5. Conclusions 

The massive star-forming region NGC7538-IRS 1 has been 
observed at 22-GHz in full polarization spectral mode with 
the VLBA and at 6.7-GHz with the EVN and MERLIN to 
detect linear and circular polarization emission from H2O and 
CH3OH masers, respectively. We detected 17 H2O masers and 
49 CH3OH masers at high angular resolution. We have mea- 
sured Zeeman-splitting for three CH3OH masers ranging from 
-2.7 ms"' to -1-2.7 ms"'. No significant magnetic field strength 
has been measured from the H2O masers. Furthermore, we have 
also shown that the masers of NGC7538-IRS 1 are all consistent 
with a torus-outflow scenario. Here the CH3OH masers are 
tracing the interface between the infall and the large-scale torus, 
and the H2O masers are related to the blue-shifted part of the 
outflow. The H2O masers of the southern group are thought 



to be associated with another source. The magnetic field is 
situated on the two surfaces of the torus with a counterclockwise 
direction on the top surface. 
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